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Abstract 

We discuss dual time evolution scenarios which, albeit running according to 
the same real time clock, in each considered case may be mapped among each 
other by means of a suitable analytic continuation in time procedure. This dy- 
namical duality is a generic feature of diffusion-type processes. Technically that 
involves a familiar transformation from a non-Hermitian Fokker-Planck operator 
to the Hermitian operator (e.g. Schrodinger Hamiltonian) , whose negative is 
known to generate a dynamical semigroup. Under suitable restrictions upon the 
generator, the semigroup admits an analytic continuation in time and ultimately 
yields dual motions. We analyze an extension of the duality concept to Levy 
flights, free and with an external forcing, while presuming that the correspond- 
ing evolution rule (fractional dynamical semigroup) is a dual counterpart of the 
quantum motion (fractional unitary dynamics). 

PACS numbers: 02.50.Ey, 05.20.-y, 05.40.Jc 

1 Brownian motion inspirations 

1.1 Diffusion- type processes and dynamical semigroups 

The Langevin equation for a one- dimensional stochastic diffusion process in an external 
conservative force field F = — (VV): x = F(x) + V2Db{t), where b{t) stands for the 
normalized white noise {b(t)) = 0, {b{t')b{t)) = 6(t — t'), gives rise to the corresponding 
Fokker-Planck equation for the probability density 'p{x,t): 

dt^ = DAp-V{Fp). (1) 

By means of a standard substitution 'p{x,t) = '^{x,t) exp[—V{x) /2D], pL], we pass to 
a generalized diffusion equation for an auxiliary function \l/(x, t): 

dt^ = DA^ - V{x)^ (2) 
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where a compatibility condition V(x) = (1/2)[(F^/2D) + VF] needs to be respected. 
This transformation assigns the role of the dynamics generator to the Hermitian (even- 
tually self-adjoint) operator —H = DA — V. 

Under suitable restrictions upon V{x), —H becomes a legitimate generator of a 
contractive dynamical semigroup exp(— ^t), t > 0. If additionally the dynamical 
semigroup is amenable to an analytic continuation in time, the contractive semigroup 
operator exp(— ift) can be related with the unitary operator exp{—iHt) via so-called 
Wick rotation t it. This duality observation underlies our forthcoming discussion 
and generalizations to Levy flights framework. 

1.2 Free propagation and its analytic continuation in time 

The standard theory of Gaussian diffusion-type processes takes the Wiener process as 
the "free noise" model, with the Laplacian as the "noise" generator. It is an element 
of folk lore that the related dissipative semigroup dynamics exp(tDA) = exp(— ti^o) 
(and thus the heat equation) can be mapped into the unitary dynamics exp{itDA) = 
exp{—itHo) (and thus the free Schrodinger equation), by means of an analytic con- 
tinuation in time procedure, [2|. A parameter D may be interpreted dimensionally as 
D = h/2m, OT D = kBT/mP (Einstein's fluctuation-dissipation statement). 

Quite often, this mapping is represented by a formal it ^ t time transformation of 
the free Schrodinger picture dynamics (one should be aware that to execute a mapping 
for concrete solutions, a proper adjustment of the time interval boundaries is necessary): 

idt^ = -D/\i} — > dtO^ = DAe^ , (3) 

where the notation 6^ for solutions of the heat equation has been adopted, to stay in 
conformity with the forthcoming more general discussion, where 6^{x,t) needs not to 
be a probability density, [2]-[l]. 

The mapping is usually exemplified in terms of integral kernels g and k as follows, 
c.f. also [5]: 

ip{x,t) = I dx' g{x — x' ,t)ip{x\0) (4) 



g{x - x', t) = k{x - x', it) = {AniDty^^^ exp[- ^^^^^^^ ] 

and 

9^{x,t) = j dx'k{x-x',t)e^{x',0) (5) 

k{x - x', t) = g{x - X, -it) = (AnDt)^/^ exp[- ^^^^ ^ ] , 

where the initial t = data need to be properly adjusted. Here, g{x — x',t) is an 
integral kernel of the unitary evolution operator: [exp[iDtA) iIj]{x,0) = ilj{x,t). The 
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heat kernel k{x — x' ,y) plays the same role with respect to the contractive semigroup 
operator: [exp{DtA) 9^]{x,0). 
The special choice of 

„2 



ip{x, 0) = (7ra2)-^/'^exp 
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0^{x,t) = ip{x, -it) = ( — ) {a^ + 2Dt)-^'^ex^ 



(6) 
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with e^{x,0) = ilj{x, 0). 

We note that p — j^'P = i/^ip* is a quantum mechanical probabihty density on R 
for all times 



p(x,t) = 



1/2 



exp 



a'^x'^ 



a" 



(9) 



The real solution 9^{x,t) of the heat equation is not a probability density p{x,t) = 
9^.{x, t) 9{x, t), unless multiplied by an appropriate real function 9(x, t) which solves the 
time adjoint heat equation (that becomes an ill-posed dynamical problem if considered 
carelessly) . 

Case 1: Since p{x,t) = [2n{a'^ + 2Dt)]-'^/'^ expl-x"^ /2{a'^ + 2Dt)] actually is an 
example of the free Brownian motion probability density for all t > 0, we infer 



p{x,t) = {47ia^y/^94x,t) = {99,){x,t) 



(10) 



where 9{x,t) = 9 — {Ana'^Y^'^ is interpreted as a trivial (constant) solution of the time 
adjoint heat equation dt9 — —DA9. We stress that 9* — (47rQ;^)'~^/^p ~ p. This, 
looking redundant observation, will prove quite useful in a more general framework to 
be introduced in below. 

Case 2: A complex conjugate t/j*{x,t) — ip{x,—t) of 'ip{x,t), Eq. (7), solves the 
time-adjoint Schrodinger equation idtip* — DAip*. Hence a time-symmetric approach 
to the analytic continuation in time might look more compelling. Indeed 



9{x,t) ^ip*{x,it) 



a 
n 



2\ 1/4 



[a^ - 2Dt)-^/^ exp 



X 



2{pP- - 2Dt) 



(11) 



is a legitimate solution of the time-adjoint heat equation dt9 — —DA9 as long as 
t e [-T/2, +T/2] where T = a'^/D. 

In the present case, both time adjoint equations set well defined Cauchy problems 
(at least in the just defined time interval). The subtle point is that the would be 
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"initial" data for the backward in time evolution, in fact need to be the terminal data, 
given at the end-point T/2 of the considered time- interval. 

The only propagation tool, we have in hands, is the heat kernel (3): k{x — — >■ 
t - t') with t > t'. There holds e^{x, t) = J k{x - x',t- t') 9^{x', t') dx' and e{x', t') = 
J e{x, t) k{x -x',t- t') dx for any t' < t e [-T/2, +T/2]. 

The original quantum mechanical probability density p = lipl"^ = ipip*, Eq. (7), is 
mapped into a Brownian bridge (pinned Brownian motion) probability density: 



p{x,±it) = p{x,t) = {99^){x,t) 



1/2 



exp 



a^x"^ 



(12) 



The price paid for the time-symmetric appearance of this formula is a limitation of the 
admissible time span to a finite time-interval of length T = / D. 

Case 3: To make a direct comparison of Case 2 with the previous Case 1, let us 
confine the time interval of Case 2 to [0, +T/2]. Now, a conditional Brownian motion 
connects p{x,0) = p(x, 0) = (a'^n)-^/^ exp{-x'^ /a"^) with p(x,t +T/2) of Eq. (10). 
Because of T = a'^/D, as t — T/2, instead of a regular function we arrive at the linear 
functional (generalized function), here represented by the Dirac delta 6{x). Note that 
6{x — x') is a standard initial t = value of the heat kernel k{x — x', t). 

This behavior is faithfully reproduced by the time evolution of 9^{x,t) and 6{x,t) 
that compose p{x,t) = {6^6){x^t) for t G [0,T/2]. The initial value of 6'*(x,0) = 
'?/'(x, 0), Eq. (6), is propagated forward in accordance with Eq. (8) to 9^:{x,T/2) = 
{Ana^y^/^expi-xyAa^). 

In parallel, 9{x,t) of (11) interpolates backwards between 9{x,T/2) = (Aira'^Y^^ S{x) 
and 9{x, 0) = 9^{x, 0). We have here employed an identity 6{ax) = {l/\a\)6{x). Because 
of fix)S{x) = /(0)5(x), we arrive at p(x,T/2) = {9^ 9){x,T/2) = 6{x). 



1.3 Schrodinger's boundary data problem 

The above discussion provides particular solutions to so-called Schrodinger boundary 
data problem, under an assumption that a Markov stochastic process which interpolates 
between two a priori given probability densities p{x, 0) and p{x, T/2) can be modeled by 
means of the Gauss probability law (e.g. in terms of the heat kernel). That incorporates 
the free Brownian motion (Wiener process) and all its conditional variants, Brownian 
bridges being included, [3l H] and [6]- [8], c.f. also [2J. 

For our purposes the relevant information is that, if the interpolating process is 
to display the Markov property, then it has to be specified by the joint probability 
measure {A and B are Borel sets in R): 

m{A, B) = dx dym{x,y) (13) 
J A Jb 
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where fj^m{x,y)dy = p{x, 0), and f^m{x,y)dx = p{y,T/2). From the start, we assign 
densities to all measures to be dealt with, and we assume the functional form of the 
density m{x,y) 

mix, y) = f{x)k{x, 0, T/2)g{y) (14) 

to involve two unknown functions f{x) and g{y) which are of the same sign and nonzero, 
while k{x, s, t) is any bounded strictly positive (dynamical semigroup) kernel defined 
for all times < s < t < T/2. For each concrete choice of the kernel, the above 
integral equations are known to determine functions /(x), g{y) uniquely (up to constant 
factors). 

By denoting 9^:{x,t) = J f{z)k{z,0,x,t)dz and 9{x,t) = J k{x,t, z,T/2)g{z)dz it 
follows that 

p(x, t) = 9{x, t)6^ (x, t) = p{y, s, X, t)p{y, s)dy, (15) 



p{y, s, X, t) = k{y, s, x, t)— -, 

0{y,s) 

for all < s < t < T/2. The above p{y, s,x,t) is the transition probability density of 
the pertinent Markov process that interpolates between p(x, 0) and p(x, T/2). Cases 1 
through 3 are particular examples of the above reasoning, once k{x, s, y, t) is specified 
to be the heat kernel (3) and the corresponding boundary density data are chosen. 
Clearly, e*{x,0) = f{x) while e{x,T/2) = g{x). 

We recall that in case of the free evolution, by setting 6{x, t) = 6 = const, as in Case 
1, we effectively transform an integral kernel k of the L}{R) norm-preserving semigroup 
into a transition probability density p of the Markov stochastic process. Then 9* ~ p. 

2 Free noise models: Levy flights and fractional 
(Levy) semigroups 

The Schrodinger boundary data problem is amenable to an immediate generalization 
to infinitely divisible probability laws which induce contractive semigroups (and their 
kernels) for general Gaussian and non-Gaussian noise models. They allow for various 
jump and jump-type stochastic processes instead of a diffusion process. 

A subclass of stable probability laws contains a subset that is associated in the 
literature with the concept of Levy flights. At this point let us invoke a functional 
analytic lore, where contractive semigroup operators, their generators and the pertinent 
integral kernels can be directly deduced from the Levy-Khitchine formula, compare e.g. 

m- 

Let us consider semigroup generators (Hamiltonians, up to dimensional constants) 
of the form H = F{p), where p = —iV stands for the momentum operator (up to the 
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disregarded h or 2mD factor) and for —oo < k < +00, the function F = F[k) is real 
valued, bounded from below and locally integrable. Then, 

/+00 
exp[-tF{k)]dE{k) (16) 
■00 

where t > and dE{k) is the spectral measure of p. 
Because of 

1 

{E{k)f){x) = — = / exp{zpx)f{p)dp (17) 

V^TT J -co 

where / is the Fourier transform of /, we learn that 

/+00 
eM-tF(k))dE(k)f](x) = (18) 
-00 

-j=j exp[-tF{k)] — [j exp{ipx)fip)dp]dk= (22) 



1 



+00 

^ exp{-tF{k)) exp{ikx)f{k)dk = [exp{-tF{p))f{p)]'^{x) 

27r J —00 

where the superscript V denotes the inverse Fourier transform. 
Let us set 

kt^^[exp{-tF{p)r. (19) 



Then the action of exp(— t^) can be given in terms of a convolution (i.e. by means of 
an integral kernel kt = k{x — y,t) — k{y, 0, x, t)): 

exp{-tH)f = [exp{-tF{p))f{p)Y = f * h (20) 

where 

{f*g){x):= [ g{x-z)f{z)dz. (21) 
Jr 

We shall restrict considerations only to those F{p) which give rise to positivity 
preserving semigroups: if F{p) satisfies the celebrated Lcvy-Khintchine formula, then 
kt is a positive measure for all t > 0. The most general case refers to a combined 
contribution from three types of processes: deterministic, Gaussian, and the jump- 
type process. 

We recall that a characteristic function of a random variable X completely deter- 
mines a probability distribution of that variable. If this distribution admits a density we 
can write E[exp{ipX)] = J^p{x) exp{ipx)dx which, for infinitely divisible probability 
laws, gives rise to the Levy- Khint chine formula 

/+00 
ipy 
[exp{ipy) — 1 — -]i'{dy)} (22) 
-00 1 ~l~ y 
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where I'idy) stands for the so-caUed Levy measure. In terms of Markov stochastic 
processes aU that amounts to a decomposition of Xt into 

Xt^at + aBt + Jt + Mt (23) 

where Bf stands for the free Brownian motion (Wiener process), Jt is a Poisson process 
while Mt is a general jump-type process (more technically, martingale with jumps). 

By disregarding the deterministic and jump-type contributions in the above, we are 
left with the Wiener process X( = aBt. For a Gaussian p(a;) = (27r(j^)~^/^ exp(— a;^/2(j^) 
we directly evaluate E[exp{ipx)] = cxp(— a^p^/2). 

Let us set = 2Dt. We get E [cxp{ipXt)] = cxp(— tDp^) and subsequently, by 
employing p ^ p = — iV, we arrive at the contractive semigroup operator cxp(tDA) 
where the one-dimensional Laplacian A = d? /dx^ has been introduced. That amounts 
to choosing a special version of the previously introduced Hamiltonian H = F{p) = 
Dp^. Note that we can get read of the constant D by rescaling the time parameter in 
the above. 

Presently, we shall concentrate on the integral part of the Lcvy-Khintchine for- 
mula, which is responsible for arbitrary stochastic jump features. By disregarding the 
deterministic and Brownian motion entries we arrive at: 

F(P) = - r^[exp(ipy) - 1 - .r^Hdy) (24) 

where i'{dy) stands for the appropriate Levy measure. The corresponding non-Gaussian 
Markov process is characterized by 

E[eMipXt)] = exp[-iF(p)] (25) 

with F{p), (22). Accordingly, the contractive semigroup generator may be defined as 
follows: F{p) = H. 

For concreteness we can mention some explicit examples of non-Gaussian Markov 
semigroup generators. F{p) = jIpI^ where < 2 and 7 > stands for the intensity 
parameter of the Levy process, upon p ^ p = — iV gives rise to a pseudo-differential 
operator H — ^N^l"^ often named the fractional Hamiltonian. Note that, by construc- 
tion, it is a positive operator (quite alike —DA). 

The corresponding jump-type dynamics is interpreted in terms of Levy flights. In 
particular 

F{V) = l\p\ F{p) = 7|V| = 7(|A|)^/' (26) 

refers to the Cauchy process. 

Since we know that the probability density of the free Brownian motion is a solution 
of the Fokker-Planck (here, simply - heat) equation 

dtp = DAp (27) 
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it is instructive to set in comparison the pseudo-differential Fokker-Planck equation 

wliicli corresponds to tlie fractional Hamiltonian and the fractional semigroup exp{—tH) = 
exp(-7|A|^/2) 

dtp = -^{Al^/^p . (28) 

As mentioned in the discussion of Case 1, instead of p in the above we can insert 6^, ~ p, 
while remembering that 6 = const. 

3 Free fractional Schrodinger equation 

Fractional Hamiltonians H = 7|A|'^/^ with p < 2 and 7 > are self-adjoint operators 
in suitable L'^{R) domains. They are also positive operators, so that the respective 
fractional semigroups are holomorphic, i. e. we can replace the time parameter t by a 
complex one a = t + is, t > so that 

exp{-aH) = / exp{-aF{k)) dE{k) . (29) 
Jr 

Its action is defined by 

[exp(-a^)]/ = [(/exp(-aF)]^ = f*K. (30) 

Here, the integral kernel reads k„ = ■^^[exp(— crF)]^. Since H is selfadjoint, the 
limit t I leaves us with the unitary group exp{—isH), acting in the same way: 
[exp{—isH)]f = [fexp{—isF)\^, except that now his := -^[exp{—isF)]'^ in general is 
not a probability measure. 

In view of unitarity, the unit ball in is an invariant of the dynamics. Hence 
probability densities, in a standard form p = ip* ip can be associated with solutions of 
the free fractional (pseudodiferential) Schrodinger equations: 

iatV(a;,t) = 7|A|'^/2^(x,t) (31) 

with initial data iIj{x,0). Attempts towards formulating the fractional quantum me- 
chanics can be found in Refs. [H [TTl [121 fT3j . 

All that amounts to an analytic continuation in time, in close affinity with the 
Gaussian pattern (1): 

idtij = 7|A|^/2^ i — > dtO* = -7|A|^/2^* (32) 

We assume that 6* ^J) and thence the corresponding 6 = const. 

Stable stochastic processes and their quantum counterparts are plagued by a com- 
mon disease: it is extremely hard, if possible at all, to produce insightful analytic 
solutions. To get a flavor of intricacies to be faced and the level of technical difficulties. 
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uoidy) := limtio[-k{0,0,y,t)]dy = — . (35) 



we shall reproduce some observations in regard to the Cauchy dynamical semigroup 
and its unitary (quantum) partner. For convenience we scale out a parameter 7. 

For the Cauchy process, whose generator is | V|, we deal with a probabilistic classics: 

'^ = 1 ^ ^' = l {t-s)'l{x-yy ^ 

where < s < t. We have {exp[ipX (t)]) := Jj^exp{ipx)p{x,t) dx = exp[— tF(p)] = 
exp(— |p|t) and 

p{x,t) = / k{y,s,x,t)p{y,s)dy (34) 

JR 

for alH > s > 0. We recall that Urritio = S{x). 

The characteristic function of k{y, s,x,t) for y, s fixed, reads exp[ipy — \p\{t — s)], 
and the Levy measure needed to evaluate the Levy-Khintchine integral reads: 

hiO,0,y,t)]dy = ^^ 
t ny^ 

To pass to a dual Cauchy-Schrodinger dynamics, we need to perform an analytic 
continuation in time. We deal with a holomorphic fractional semigroup exp(— (Tt| V|), 
a = t + is, (27). It is clear that exj9(— 1| V|) and exp{—is\V\) have a common, identity 
operator limit as t | and s = t | 0. 

An analytic continuation of the Cauchy kernel by means of (28) gives rise to: 

It 1 1 is 

^*(^) = -^9s{x) = kUx) = -[S{x -s) + 6{x + s)] + -V- , (36) 

7C X^ 2 IT X^ — S^ 

where V indicates that a convolution of the integral kernel with any function should be 
considered as a principal value of an improper integral, [8]. This should be compared 
with an almost trivial outcome of the previous mapping (2)—* (3). Here, we employ the 
usual notation for the Dirac delta functionals, and the new time label s is a remnant 
of the limiting procedure t I in a = t + is. 

The function denoted by is/Tc{x'^ — s^) comes from the inverse Fourier transform of 
— ^^sgn{p)sin{sp). Because of 

[sgn{p)r = i\I^V{-) (37) 

V TT X 

where 'P(^) stands for the functional defined in terms of a principal value of the integral. 
Using the notation 5±s for the Dirac delta functional 5{x =F s)\ 

[sin{sp)Y = i^{5,-5^s) (38) 

we realize that 

^ = ^(5.-M*^(-) (39) 

TT — 27r X 



9 



is given in terms of the implicit convolution of two generalized functions. Obviously, a 
propagation of an initial function ipoi^) to time t > 0: 

'ijj{x,t) = / g{x — x' ,t)'ijjo{x')dx' (40) 
Jr 

gives a solution of the fractional (Cauchy) Schrodinger equation idtip = — |V|V'- 

In comparison with the Gaussian case of Section 1, one important difference must be 
emphasized. The improper integrals, which appear while evaluating various convolu- 
tions, need to be handled by means of their principal value. Therefore, a simple it ^ t 
transformation recipe no longer works on the level of integral kernels and respective ijj 
and 6* functions. 

One explicit example is provided by the incongruence of (31) and (34) with respect 
to the formal t —>■ —it mapping. Another is provided by considering specific solutions 
of pseudo-differential equations (30). 

To that end, let us consider 6^*0(3;) = (2/7r)^/^Y:^, together with 6 = (27r)~^/^. 
Then, 9 9^:{x,0) = ^^-^^^2) is an L{R) normalized Cauchy density, while 9^:o{x) itself is 
the L'^{R) normalized function. Clearly: 

e,{x,t) = [exp(-t|V|)^.o](a:) = J k{y,0,x,t)e4y,0)dy = 0^ ^2^"^(j^^)2 (^1) 
while the corresponding ip{x,t) with ipo^x) = 9^o{x) reads (for details see e.g. [8]): 

ipix, s) = [exp(-zs|V|) ipolix) = ^[M^ + - (42) 

-[{x - s)iIjo{x - s) -{x + s)iIjo{x + s)] . 

4 Dynamical duality in external potentials: frac- 
tional Schrodinger semigroups and Levy flights 

4.1 Schrodinger semigroups for Smoluchowski processes 

Considerations of Section 1, where the free quantum dynamics and free Brownian mo- 
tion were considered as dual dynamical scenarios, can be generalized to an externally 
perturbed dynamics, [2j. Namely, one knows that the Schrodinger equation for a quan- 
tum particle in an external potential V{x), and the generalized heat equation are con- 
nected by analytic continuation in time, known to take the Feynman-Kac (holomorphic 
semigroup) kernel into the Green function of the corresponding quantum mechanical 
problem. 

idt^ = -DA^ + V^i — > dtO, = DM, - VO, . (43) 
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Here V = V{x)/2mD. 

For V = V{x),x G R, bounded from below, the generator H = —2mD'^A + V 
is essentially selfadjoint on a natural dense subset of L^, and the kernel k{x,s,y,t) = 
[exp[—(t — s)H]]{x, y) of the related dynamical semigroup exp(— tff) is strictly positive. 
The quantum unitary dynamics exp{—iHt) is the an obvious result of the analytic 
continuation in time of a dynamical semigroup. 

Assumptions concerning the admissible potential may be relaxed. The necessary 
demands are that H is self-adjoint and bounded from below. Then the respective 
dynamical semigroup is holomorphic. 

The key role of an integral kernel of the dynamical semigroup operator has been 
elucidated in formulas (11)-(13), where an explicit form of a transition probability 
density of the Markov diffusion process was given. We have determined as well the 
time development of 9^:{x,t) and 9{x,t), so that p{x,t) = {60^){x,t) is a probability 
density of the pertinent process. 

If we a priori consider 6{x,t) in the functional form 6{x,t) = exp$(x,t), so that 
t) = t) exp[— $(x)], and properly define the forward drift b{x, t) = 2D'V^{x, t) 
in the pertinent Fokker-Planck equation: 

dtp = DAp-\/{bp) (44) 

we can recast a diffusion problem in terms of a pair of time adjoint generalized heat 
equations 

dt.e, = DAe, - ve, (45) 

and 

dtO = -DAB + V9 , (46) 

i. e. as the Schrodinger boundary data problem, where an interpolating stochastic 
process is uniquely determined by a continuous and positive Feynman-Kac kernel of 
the Schrodinger semigroup exp(— ti^), where H = —DA + V. 

If our departure point is the Fokker-Planck (or Langevin) equation with the a priori 
prescribed potential function $(x,t) for the forward drift b{x,t), then the backward 
equation (44) becomes an identity from which V directly follows, in terms of $ and its 
derivatives, [6l [7] : 

V{x,t) = [d,<l> + ^i^ + Vb)] (47) 

For the time-independent drift potential, which is the case for standard Smoluchowski 
diffusion processes, we get (c.f. also [1], where the a transformation of the Fokker- 
Planck equations (42) into an associated Hermitian problem (43) is described in detail): 

V(a:)= [i(|^ + V6)]. (48) 
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Notice that ^(x) is defined up to an additive constant. 

To give an example of a pedestrian reasoning based on tlie above procedure in 
case of a concrete Smolucliowski diffusion processes, let us begin from the Langevin 
equation for the one-dimensional stochastic process in the external conservative force 
field F{x) = —{W){x) (to keep in touch with the previous notations, note that $ = 

-vy. 

(1 T I 

^ = Fix) + V2DBit) (49) 
at 

where B{t) stands for the normalized white noise: {B{t)) = 0, {B{t')B{t)) = 6{t — t'). 
The corresponding Fokker-Planck equation for the probability density p{x, t) reads: 

dtp = DAp-V{Fp) (50) 

and by means of a substitution p{x,t) = 9^{x,t) exp[—V{x) /2D], [Ij, can be trans- 
formed into the generalized diffusion equation for an auxiliary function 6*{x,t): 

dtd, = DM, - Vd, (51) 

where the consistency condition (reconciling the functional form of V with this for F) 

^=\i^ + ^F]- (52) 



2 \2D 

directly comes from the time-adjoint equation 

dte = Q = -DM + Ve (53) 

with e{x) = exp[-V{x)/2D]. 

For the Ornstein-Uhlenbeck process b{x) = F{x) = —kx and accordingly 

2 2 

V(a;) = (54) 
^ ' AD 2 ^ ' 

is an explicit functional form of the potential V, present in previous formulas (41)- (44). 

4.2 Fractional semigroups and perturbed Levy flights 

External perturbations in the additive form: 

idti){x, t) = 7| A|^/2^(a;, t) + V{x)^ix, t) (55) 



were considered in the framework of fractional quantum mechanics, [Tl]-[T3], c.f. also 
[HI [9]. With the dual dynamics concept in mind, Eq. (30), we expect that an anlytic 
continuation in time (if admitted) takes us from the fractional Schrodinger equation to 
the fractional analog of the generalized diffusion equation: 

dtd* = -7|A|'^/2^* - Vr . (56) 
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The time- adjoint equation has the form 

dtO = ^\A\^'/^e + Ve, (57) 

We shall be particularly interested in the time-independent 6{x,t) = 6{x), an assump- 
tion affine to that involved in the passage from (44)- (46). 

Hermitian fractional problems of the form (48) and/or (49) have also been studied in 
Refs. [lll[l5l[T6]. However, the major (albeit implicit, never openly stated) assumption 
of Refs. [m [151 [16] was to consider the so-called step Levy process instead of the jump- 
type Levy process proper. 

This amounts to introducing a lower bound on the length of admissible jumps: 
arbitrarily small jumps are then excluded. That allows to by-pass a serious technical 
obstacle. Indeed, for a pseudo-differential operator 7A^/^, the action on a function from 
its domain can be greatly simplified by disregarding jumps of length not exceeding a 
fixed e > 0, see e.g. Refs. [8l[9]: 

l\A\^/'f){x) = - I [f(x + y)-f{x)-y^P$]u{dy) (58) 
Jr 1 + 1/ 

7|A|^/V)(x) = - / [f{x + y)-f{xMdy). 

Compare e.g. Eq. (2) in [15J and Eq. (6) in [TB]. Note that these Authors have skipped 
the minus sign that must appear on the right-hand-side of both formulas (50). 

As a side comment, let us point out that the principal integral value issues of Section 
3 would not arise in our previous discussion of Cauchy flights and their generators, if 
arbitrarily small jumps were eliminated from the start. Nonetheless, if the e | limit 
is under control, the step process can be considered as a meaningful approximation 
of the fully-fiedged (perturbed) jump-type Levy process. This approximation problem 
has been investigated in detail, in the construction of the perturbed Cauchy process, 
governed by the Hermitian dynamical problem (53), with the input (55), under suitable 
restrictions on the behavior of V, |9] . 

Let us come back to time-adjoint fractional equations (54) and (55). We have 
p(x,t) = {66*){x,t) and employ the trial ansatz of Section 4.2: 

e{x,t) = 9{x) =exp[^{x)] (59) 

9*{x, t) = p(x, t) exp[-$(x)] . 
Accordingly (55) implies, compare e.g. [H] for an independent argument: 

V = —7 exp (_$)|A|^/2exp($) (60) 
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to be compared with Eq. (8) in Ref. [IS]. In view of (54) we have 



d{p = edtO* = -7exp(0)[| Al'^/^ exp(-$)p] + Vp = -Vj . 



(61) 



Langevin-style description of perturbed Levy flights (deterministic component plus 
the Levy noise contribution) are known, [TTl UHl HH], to generate fractional Fokker- 
Planck equations of the form 



where F = —W = V$, we face problems which are left unsettled at the present stage 
of our investigation: 

(i) May the stochastic processes driving (58) and/or (59) coincide under any cir- 
cumstances, or basically not at all ? 

(ii) Give an insightful/useful definition of the probability current j{x,t) in both 
considered cases, while remembering that for fractional derivatives the composition 
rule for consecutive (Riesz) derivatives typically breaks down. 

Both problems (i) and (ii) have have an immediate resolution in case of diffusion- 
type processes, where by departing from the Langevin equation one infers Fokker- 
Planck and continuity equations. In turn, these equations can be alternatively derived 
by means of the Schrodinger boundary data problem,provided its integral kernel stems 
from thee Schrodinger semigroup, both in the free and perturbed cases. The stochastic 
diffusion process (corresponding to that associated with the Langevin equation) is then 
reconstructed as well. Thence, the Schrodinger loop gets closed. 

While passing to Levy processes, we have demonstrated that, with suitable reserva- 
tions, this Schrodinger "loop" can be completed in case of free Levy flights. However, 
the "loop" remains incomplete (neither definitely proved or disproved) for perturbed 
Levy flights. 

At this point we should mention clear indications [H] that, once discussing Levy 
flights, we actually encounter two different classes of processes with incompatible dy- 
namical properties. One class is related to the Langevin equation, another - termed 
topological - relies on the "potential landscape" provided by the effective potential 
V{x). An extended discussion of the latter problem has been postponed to the forth- 
coming paper, c.f. [23] . 

Acknowledgement: Partial support from the Laboratory for Physical Founda- 
tions of Information Processing is gratefully acknowledged. 



dtp = -V(Fp) - 7| Al^/^p = -Vj 



(62) 



14 



References 

[1] H. Risken, The Fokker-Planck equation, Sringer-Verlag, Berlin, 1989 
[2] P. Garbaczewski, Phys. Rev. E 78, 031101, (2008) 
[3] J-C. Zambrini, J. Math. Phys. 27, 2307. (1986) 
[4] J. C. Zambrini, Phys. Rev. A 35, 3631, (1987) 
[5] M. S. Wang, Phys. Rev. A 37, 1036, (1988) 

[6] Ph. Blanchard and P. Garbaczewski, Phys. Rev. E 49, 3815, (1994) 

[7] P. Garbaczewski and R. Olkiewicz, J. Math. Phys. 37, 732, (1996) 

[8] P. Garbaczewski, J. R. Klauder and R. Olkiewicz, Phys. Rev. E 51, 4114, (1995) 

[9] P. Garbaczewski and R. Olkiewicz, J. Math. Phys. 40, 1057, (1999) 

[10] P. Garbaczewski and R. Olkiewicz, J. Math. Phys. 41, 6843, (2000) 

[11] N. Laskin, Phys. Rev. E62,3135, (2000) 

[12] N. Laskin, Phys. Rev. E66, 056108, (2002) 

[13] N. Cufaro Petroni and M. Pusterla, Physica A 388, 824, (2009) 

[14] D. Brockmann and I. Sokolov, Chem. Phys. 284, 409, (2002) 

[15] D. Brockmann and T. Geisel, Phys. Rev. Lett. 90, 170601, (2003) 

[16] D. Brockmann and T. Geisel, Phys. Rev. Lett. 91, 048303, (2003) 

[17] S. Jespersen, R. Metzler and H. C. Fogedby, Phys. Rev. E 59, 2736, (1999) 

[18] P. D. Ditlevsen, Phys. Rev. E 60, 172, (1999) 

[19] V. V. Janovsky et al., Physica A 282, 13, (2000) 

[20] A. Chechkin et al., Chem. Phys. 284, 233, (2002) 

[21] A. A. Dubkov and B. Spagnolo, Acta Phys. Pol. A 38, 1745, (2007) 

[22] A. A. Dubkov, B. Spagnolo and V. V. Uchaikin, Int. J. Bifurcations and Chaos, 
18, 2549, (2008) 

[23] P. Garbaczewski, Levy flights and Levy-Schrodinger semigroups, larXiv:0902. 35361 



15 



